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W-PARTICLE DISTRIBUTION IN ELECTROWEAK 
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Valckenierstraat 65, 1018 XE Amsterdam, the Netherlands 

Results are presented of a numerical study of the distribution of W-bosons gener- 
ated in a tachyonic electroweak prc-heating transition. 

1. W-distribution in tachyonic electroweak baryogenesis 

The electroweak transition may have taken place shortly after inflation 
ending at a low energy scale, and this possibility has been used to suggest 
alternative scenarios for baryogenesis. We have studied such a transition 
in the SU(2)-Higgs model with effective CP-violation . The transition 
is assumed to have taken place at zero temperature. In the quenching 
approximation, it is induced simply by flipping the sign of the quadratic 
term in the Higgs potential, /x 2 </>t(/> — > — /i 2 0t</>. This causes the Higgs field 
to go through a spinodal instability in which its particle numbers initially 
grow exponentially fast, leading to classical behavior 2 > 3 ' 4 > 5 ' 6 > 7 . The gauge 
fields react strongly, and it is interesting to see the emergence of effective 
W-particles, their energy spectrum u>k, and their distribution function rife. 

The numerical simulation is carried out in the temporal gauge Aq = 0. 
To define LOk and rik, we transform to the Coulomb gauge djAj — (which is 
a smooth gauge in which it makes sense to Fourier transform to momentum 
space) , and maximally coarse-grain the equal-time correlators of the gauge 
field and its canonical conjugate (Ej) over the (periodic) volume L 3 , 
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where we also averaged over the three isospin directions. After Fourier 
transformation, averaging over the transverse modes, 

\Saj d 3 xe-^C AA (x)=C AA , 

2 (j**J — ^ikjj J d?x e * k x C EE (x) = C EE , 

and averaging over directions, k — > k = |k|, the time-dependent particle 
spectrum and distribution function is obtained by solving 

C AA = n k /ia k C EE = n k u k 

for n k and ui k . Note that n k appears and not n k + 1/2, since we are using 
the classical approximation. 




We present here results obtained from only one gauge field configu- 
ration, i.e. the average over initial conditions (• • •) is not carried out. If 
classical behavior applies, then one field configuration corresponds to one 
realized universe. Of course, statistical observables such as n k and uj k will 
suffer from 'cosmic variance', which may be severe in a small volume. The 
parameters of the simulation on a lattice of 60 3 sites are given by 

couplings: g 2 = 4/9, I = 8X/g 2 = m^/m^ 
physical size L = 21m», lattice spacing a = 0.35mj 1 1 

Figure 1 shows n k at early times after the quench at t = 0. We see 
an exponential growth of the low momentum modes until they saturate 
around t ~ 8 mj^ , after which the tail quickly acquires a different slope on 
the log-plot. Figure 2 shows uj k for t < 20 . It suffers greatly from 
fluctuations ('cosmic variance' in our small 'universe'), and we averaged 
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Figure 2. a 2 ^ 2 , versus a?k 2 at early times. The particle interpretation starts to make 
sense from tmjj 18 onwards. 
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Figure 3. i 2 ^| versus a 2 fc 2 at later times. The fitted slopes are 1 within a few percent. 
The fitted intercepts correspond to effective masses m^s/mjj as 0.86, 0.88, 0.95, 0.85 
respectively for tm H = 100, 200, 400, 1000. 

the data in bins of size A|k| 0.14 m^. Only from t = 18 mjj 1 onwards 
a sensible u>k emerges, as can also be seen at later times in Figure 3. The 
spectrum is close to the simple form u>\ = m 2 s + k 2 (on the lattice we use 
k 2 — > a~ 2 J2j [2 — 2 cos(afcj)]), with m e g ~ rnw, the zero-temperature input- 
value. At later times 20 < tm,H < 1000 the distribution Uk does not change 
much, as can be seen in Figure 4. At k = 0, the particle number stays 
roughly constant, 6 < no < 7, while the tail of the distribution appears to 
rise slowly. There is little power beyond acuk = amu = 0.35 (top plot). 

The data do not yet show clear signs of classical equilibration, nk — * 
T/ujk- One expects a plateau to develop in rikWk, but the 'plateau' in the 
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Figure 4. nj. (top), n^LD^a (bottom left) and ln(l + l/n^) (bottom right) versus (Jfcd, 
at times tm H = 20, 100, 400, 1000. 

lower-left plot does not look very convincing yet as it involves only two 
fc-values. The lower-right plot in Figure 4 shows ln(l + n^ 1 ), which for a 
Bose-Einstein distribution would have the linear form (uj^ — jj,) /T . The data 
are roughly compatible with this behavior, the fits in Figure 5 indicate an 
effective temperature that slowly increases as the higher momentum modes 
get more occupied, T/m H = 0.44 — * 0.56, with a hint of a decreasing 
effective chemical potential fi « 0.8 —* 0.7, as t increases from 20 to fOOO 
m H \ 

2. Conclusion 

After t w 20 m Tf, the particles produced by the instability settle into an 
effective temperature of about 0.5 mjy. Subsequent evolution is slow in the 
SU(2)-Higgs model, in the classical approximations used here: the modes 
are massive and the effective interactions appear to be weak at the BE- 
temperature T w 0.5 mn- The Higgs particles have roughly the same effec- 
tive temperature. 
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Figure 5. Linear fits ln(l + 1/rtfc) ~ (ojj. — /u)/T resulting in an effective temperature 
T/m H = 0.44, 0.46, 0.49, 0.56, and effective chemical potential u/m H = 0.8, 0.9, 0.9, 

0. 7. respectively for tm H = 20, 100, 400 and 1000. 

The initial conditions for this simulation do not put power into the high- 
momentum modes near the cutoff ' . There is furthermore little draining of 
power from the low- to high- momentum modes, so Ray leigh- Jeans effects 
are negligible. Lattice artefacts are also under control, since the distri- 
bution function drops exponentially fast and there is little power in the 
gauge-field modes beyond ak = 0.5. In particular, it makes sense to use a 
simple lattice version 1 of an effective CP-violating term kTi <j)' (j>F tll ;F fJ " 1 ' in 
the lagrangian. 
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